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PIIEFACE 

The  object  of  this  study  is  to  determine  the  coefficient  of  eddy 
conduction,  its  variation  -rrith  height,  and  its  diurnal  variation  in  the 
surface  layer* 

This  paper  vtq.3   prepared  at  the  U,  S,  ITaval  Postgraduate  School, 
Ifonterey,  California  during  the  period  December  194?  to  Liiy  1950« 

The  author  expresses  his  appreciation  and  thanks  to  Associate 
Professor  Frar±  L«  llartin,  7/ho  gave  valuable  advice  and  guidance  in 
the  development  of  this  paper* 


(ii) 


TABLE  OF  COIjTMTS 

Page 

CmTIPICATE  OF  APPROVAL  i 

FRBFACE  ii 

TABLr]  OF  COITTEITTS  iii 

LIST  OF  ILLUSTRATIOIIS  iv 

TABLE  OF  SyUBOLS  AlTD  ABBREVIATIOIIS  v 
CHAPTER 

I.     UTTRODUGTIOK  1 

II.     TKE  IJATURE  OF  THE  EtATA  4 

1,  The  Data  4 

2,  Selection  of  Data  for  This  Study  4 

III.  THEORETICAL  AtTD  EIvlPIllICAI.  CONCEPTS  d 

!•  The  'iTater  Vapor  Radiati-ve  Transfer  Problem  6 

2»  The  Carbon  Dioxide  Radiative  Transfer  Problem  7 

5»  Heat  Transfer  by  I'olecular  Conductivity  11 

4«  Hate  of  Change  of  Temperature  hy  Eddy  Conductivity  12 

5.  The  Coefficient  of  Eddy  Conduction  12 

IV.  RATES  OF  CHMTGE  OF  TEI'I'ERATIJRE  l6 

1*  Ptato  of  Change  of  Temperature  Due  to  Yfater  Vapor  l6 
2.  Rate  of  Change  of  Temperature  Due  to  Carbon  Dioxide    22 

5.  Resultant  Radiative  Rate  of  TemperatiAre  Change  25 
4.  Rate  of  Change  of  Temperature  Due  to 

Ilolccular  Conduction  27 
5»  Rate  of  Change  of  Temperature  Due  to 

Eddy  Conduction  27 

V.  RESULTS  AITD  COITCLUSIOIIS  2? 

1»  Llethod  of  Computation  of  the  Coefficient  of 

Eddy  Conduction  29 

2.  Discussion  of  Results  41 

BIBLIOGRAPHY  47 


(iii) 


LIST  OF  ILLUSTRATIOITS 


Page 


Tables 

1.       Heat  Transfer  Due  to  V«Uter  Vapor,  cal   (5  hr)"^  20 

2*       Three-hourly  Rate  of  Cliaaige  of  lemiDorature  Due  to 

■Jater  Vapor  21 

5«       Heat  Transfer  Due  to  Carbon  Dioxide,   cal   (5  hr)  25 

4«   Tliree-hourly  Rate  of  Temperature-Change  Due  to 

Carbon  Dioxide  24 

5»   Iialf-hovirly  Rate  of  Temperature-change:  Y/ater  Vapor 

Carbon  Dioxide  26 

6«   Half -hourly  Rate  of  Temperature-change  Due  to 

Eddy  Conduction  2? 

7«   Half -hourly  Temperature -change  Due  to  Eddy  Conduction         28 

8»   Computation  Values  for  Determination  of  K  at 

Reference  Level  yi 

9»   The  Coefficient  of  Eddy  Conduction,  K,  (cm^  sec""^)  34 

10.  Values  of  Lapse  Rates  at  0200  and  IJOO  LCT  42 

11.  Slopes  of  Log  K  against  Log  Z  for  the  Lines  of  Plate  I        44 

Figxires 

1.  Difference  of  ITet  Fluxes  for  a  Slab  Below  288  Feet 

(T^-pical  ITight-time  Conditions)  17 

2.  Difference  of  ITet  Flxixes  for  a  Slab  Below  288  Feet 

(T-'pical  Day-tiiue  Conditions)  18 

5«  Adiabatic  Velocity  Profile  (A,1I,  -  P.II.  Average)  yi 

4«  Diurnal  Variation  of  K  at  .25  Feet  56 

5»  Diurnal  Variation  of  K  at  6  Feet  57 

6«  Diurnal  Variation  of  E  at  JO  Feet  58 

7f   Diurnal  Variation  of  K  at  60  Feet  39 

8*   Diurnal  Variation  of  a  at  125  Feet  40 

PLATE  (inside  back  cover) 


Variation  of  K  with  Height 


(iv) 


TABLE  OF  Syi30LS  AiHD  ABEREVIATIOHS 

A    percentage  energy  absorbed  by  a  thickness,  c,  of  carbon  dioxide 

c    the  equiralent  thickness  of  carbon  dioxide  in  cms 

c_   specific  heat  of  air  in  gn  cal 

E    flux  in  the  band  15»3  -  17 •!  microns  entering  the  second 
boundary  of  a  slab 

E    flux  in  the  band  13 .^  -  17 •!  microns  entering  the  first 
boundary  of  a  slab 

E.  black  body  flux  in  the  band  Ip*?   -  17 •!  microns  for  the  mean 
temperature  of  the  slab 

F.  net  upiTard  flux  at  bottom  of  the  slab 
F    net  upv^-ard  flux  at  top  of  the  slab 

■^m  upvra.rd  flux  of  hoat  due  to  molecular  conduction 

g  acceleration  of  gravity 

K  coefficient  of  eddy  conduction 

k  Von  Karman  constant  (k  s  •45) 

1  mixing  lengiih  in  a  non-isotropic  atmosphere 

1  mixing  length  in  an  isotropic  atmosphere 

m  slope  of  the  curve  log  K  against  log  Z 

n,x  coefficients  given  by  Callendar  for  computing  monochromatic 
absorptivity  in  the  carbon  dioxide  band  15 •O  -  17 #0  microns 

p  atmospheric  pressure  in  millibars 

p  a"bmo spheric  pressure  in  atmospheres 

p.  mean  pressure  of  the  i-th  slab 

q.  mean  specific  humidity  of  the  i-th  slab 

R ,  gas  constant  per  gram  of  dry  air 

T  absolute  temperature 

T  mean  absolute  temperature 


(v) 


u  horizontal  vfind  speed 

u  horizontal  wind  speed  in  an  isotropic  atmosphere 

Z  heifjht  measured  from  the  base  of  the  tovrer 

c^  coefficient  of  molecular  conduction 

y  the  existing  lapse  rate 

If  dry  adiahatic  lapse  rate 
d 

0  potential  teriperature 

TT  constant  stress  in  the  surface  layer 

XT  constant  stress  in  the  surface  layer  of  oji  isotropic  atmosphere 


(vi) 


I,  IKTRODUCTIOIT 

During  I948-I949  several  articles  have  "been  published  by  the 
University  of  Te;^as  giving  half -hourly  observations  of  both  tempera- 
ture and  vdnd  at  various  levels  from  the  ground  up  to  and  including 
288  feet  for  periods  of  approximately  one  day#  This  suggested  the 
problen  of  detemining  the  coefficient  of  eddy  conductivity  using  the 
observed  lapse  rate  data;  and,  the  coefficient  of  viscosity  by  using 
the  observed  -find  profile  data.  These  two  coefficients  were  found,  in 
1536*  ^y  Sverdrup   LllJ  "to  ^^  equal  in  tho  atmosphere,  but  other  eddy 
diffusion  coefficients  -were  found  by  hin  to  be  unequal.  On  the  other 
hand.  Potters  sen  and  Strinbank  I.10J  ,  in  154?  j  suggest  that  the  foriTier 
coefficient  is  more  nearly  related  to  the  latter  by  the  factor  .65  in 
tho  free  atmosphere.  The  question  appears  to  require  another  critical 
excmination;  and,  vrith  the  availability  of  new  and  fairly  accurate  data, 
it  seemed  to  be  a  problem  well  worthy  of  study.  In  actuality,  the  problem 
of  determining  the  diurnal  variation  of  the  coefficient  of  eddy  conductivity, 
K,  near  the  ground  proved  to  be  so  complex  that  this  vreis  all  that  time  per- 
mitted in  this  dissertation.  It  vrill  be  necessary  to  ans-^rer  the  question 
of  the  equality  of  the  ^ao   coefficients  at  some  later  date,  ilov/ever,  tho 
determination  of  values  of  the  coefficient  of  eddy  conductivity  and  its 
diurnal  variation,  determined  purely  from  heat  transfer  considerations, 
is  of  some  practical  interest  as  a  check  on  previously  determined  values. 
For  example,  Stev/art  ^ij  indicates  a  possible  range  of  10^  cm^  sec"l  to 
10^  cm  sec"!  for  K  within  the  surface  lajrer,  depending  upon  the  stability 
of  the  layer -» 
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This  study  gives  explicit  values  of  K,  vrith  respect  to  local  civil 
time,  corresponding  to  a  particular  synoptic  situation,  tliat  for  the 
period  27-28  September  1948  in  the  vicinity  of  I.lanor,  Te:cas«  Another 
important  result  of  this  sttidy  proved  to  be  the  determination  of  rates 
of  radiative  heating  and  coolinr  of  the  layers  of  air  near  the  ground  at 
various  times  covering  a  full  day.  This  yielded  some  surprisingly  large 
values.  These  results  could  "be  used  to  solve  a  problem  proposed  by 
Brunt  L2j  concerning  the  magnitude  of  the  coefficient,  K^,  of  radiative 
diffusivity.  Brunt  originally  computed  K.-^   to  be  about  650  cm^  soc"l  but 
in  1940,  using  the  new  absorption  coefficients  of  Dennison,  Ginsberg  and 
Weber  l4J  he  concluded  that  this  value  of  Kj^  probably  should  bo  increased 
thirty-fold*  The  Tra,tor  vapor  radiative  transfer  problem  of  this  thesis 
used  the  Dlsasser  Atmospheric  Radiation  Chart  which,  in  turn,  is  based 
primarily  on  the  new  absorption  coefficients  mentioned  above.  Actually 
K^  \ms  not  computed  in  this  study  but  the  results  of  this  study  form  a 
basis  which  vrould  permit  the  calculation  fairly  simply. 

It  ims  found  also  that  carbon  dioxide  plaj;^  a  stirprisingly  large  part 
in  the  radiative  heat  transfer  processes  near  the  ground  where  rather  large 
lapse  rates  are  fovmd.  In  the  preparation  of  his  radiation  chart,  Slsasser 
made  the  assumption  that  the  carbon  dioxide  contribution  to  the  not  flux 
could  be  neglected  in  comparison  to  that  due  to  vreiter  vapor.  It  turned 
out  in  this  study  that  the  radiative  contribution  to  the  rate  of  clionge 
of  temperature  due  to  carbon  dioxide  ivas  in  the  neighborhood  l/p  to  I/4 
that  due  to  water  vapor^  Elsasser*s  remark,  concerning  the  effect  of  carbon 
dioxide,  "that  even  thin  layers  radiate  as  black  bodies  within  these  bands" 
was  shovm  to  be  definitely  invalid.  In  fact,  using  a  formula  for  carbon 
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dioxide  absorptivity  due  to  Callendar  [^J  i  '^h©  conputed  absorptivities 
ranged  anywhere  from  2^  to  20/i  for  layers  near  the  earth's  surface,  de- 
pending upon  their  thickness • 

Surprisingly  large  values  of  turbulent  temperature-change  per  half 
hour  vrere  obtained,  in  this  study,  at  three  inches  above  the  ground*  For 
e:icample,  there  vreis  a  turbulent  rate  of  cooling  of  lO^C  per  half -hour  at 
1300  local  time*  It  ivas  found,  however,  that  oven  such  large  values  irere 
negligible  insofar  as  affecting  the  magnitude  of  the  coefficient  of  eddy 
conductivity  in  the  surface  layer.  In  other  words,  -srhenever  the  lapse 
rates  wore  significantly  different  from  dry  adiabatio,  it  vra,s  foxmd  that 
the  turbulent  temperature-change  term  of  the  partial  differential  equation 
of  eddy  conduction  vra-s  negligible  in  comparison  to  the  turbulent  heat  flux 
through  either  boundary  of  a  slab  of  air»  Thus  the  surface  layer  extended 
at  least  to  80  feet  throughout  the  period  of  this  investigation  and  during 
the  night  extended  to  about  200  feet* 

Other  results  deduced  in  this  study  are  the  following:  The  coefficient 
of  eddy  conductivity  in  the  surface  la;;,'©r  appears  to  follow  a  povrer  lav/  of 
the  form  K  =  11-^   Z^,  with  m  ^1  the  exact  value  of  m  depending  in  some  way 
upon  tho  stability*  Ubreover  in  this  layer,  K  takes  on  its  maximum  value 
at  about  I50O  LCT  and  its  minimum  value  at  about  I9OO  LCT. 
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II.     THE  HATDRE  OF  TliE  DATA 

!♦       The  Data. 

Data  utilized  in  this  research  vms  taken  from  Heport  ilo.  29  com- 
piled by  the  Electrical  Engii\eering  Research  Laboratory  of  the  University 
of  Texas  [_6j  •  This  report  contains  recorded  measurements  of  temperature 
and  vapor  pressure  for  selected  levels  from  the  groxind  up  to  and  including 
288  feet.  The  soundings  are  recorded  at  loalf -hourly  intervals  tliroughout 
the  24  hour  period  chosen  to  represent  a  particular  syaoptic  situation. 
Ileasurements  of  vfind  speed  and  direction  are  recorded  for  selected  levels 
up  to  507  feet.  The  sounding  site  is  on  level  ground  under  cultivation. 

The  authors  of  [.6 J  give  the  follov/ing  statements  of  accuracy  of  the 
data. 

(a)  The  temperatures  are  accurate  to  "iTithin  JL°C  over  the  0-50°C  range. 

(b)  The  6-ninute  averaging  period  used  to  obtain  temperatures  is  not 
sufficient  to  remove  all  turbulent  fluctuations;  hoiTever,  they  do  not  deviate 
more  than  *.2°C  from  the  30-minute  mean  temperature.  This  turbulent  fluctu- 
ation of  the  6-mintite  mean  as  compared  vrith  the  50-minute  mean  will  be  called 
"the  gust  error". 

(c)  The  Aerovanes  show  a  uniform  response  with  an  accuracy  of  t2^. 
For  a  detailed  description  of  the  data  and  the  field  installation,  the 

reader  is  referred  to  the  report  of  the  data  l^J  • 
2.   Selection  of  Data  for  this  Study. 

The  data  for  the  period  September  27-28,  I948  was  selected  for  this 
study.  A  high  pressure  cell  dominated  the  area  and  gave  conditions  of 
horizontal  homjogeniety. 
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Soundings  at  three-hourly  intervals  were  used  and  conputations  of 
heat  flux  v/ere  niade  for  the  layers  of  air  bounded  hj  the  ground,  6  and 
18  inches;  5,  6,  12,  20,  35,  55,  CO,  110,  145,  l85,  2^5  and  288  feet, 
inasDTLich  as  these  are  selected  levels  for  v/hich  the  temperature  data 
Trere  available. 
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III.  THEORETICAL  AilD  ElJPIRICilL  COIIC^^PTS 

1.   The  V;ater  Yapor  Radiative  Transfer  Problem, 

Elsasser  F 5 J  gives  a  method  for  determining  the  net  flux  duo  to 
the  imter  vapor  content  of  a  Ic-jGr   of  atmosphere  tlirough  the  use  of  his 
radiation  chart.  The  corrected  optical  depth  of  vra-ter  vapor,  u,  in  a 
column  of  air  of  unit  cross-section  is  computed  by  using  the  formula: 

U  '  -  7^  -^  fi  V/o^  Ai   -p.  (3.1) 

vrhen  only  temperature  and  moisture  data  are  available  at  given  levels, 
pressures  for  these  levels  may  be  computed  from  the  surface  pressure  by 
use  of  the  integrated  form  of  the  hydrostatic  equation: 

^  ^2  -  Rj  T  C-^-^, -^^^y  (5-2) 

For  the  evaluation  of  not  flux  into  or  out  of  a  slab  of  air,  utilizing 
areas  on  Elsasser' s  Atmospheric  Radiation  Chart,  Elsasser  recommends  the 
measurement  of  these  areas  by  means  of  a  planimeter.  The  difference  of  the 
net  fluxes  for  a  slab  is  then  given  in  gram-calories  per  three  hours  by  this 
method  and  is  converted  to  rate  of  change  of  temperature  of  the  slab  by  means 
of  the  equation: 


/r./r  _c/^Tj^ey^.  ,.,) 
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In  this  equation,  F2  ^^   '^-'^^  ^^''^  upr«ard  flux  at  the  top  of  the  slab,  while 
?-,  is  the  net  upivard  flux  at  the  bottom  of  the  slab.  For  further  details 
concerning  this  method,  reference  should  be  made  to  Elsasser*s  study.  Heat 
Transfer  by  Infrared  Radiation  in  the  Atmosphere  L^J  • 
2*   The  Carbon  Dioxide  Radiative  Transfer  Problem. 

Callendar  [^5 J  gives  formulas  for  the  absorptivity  in  terms  of  the 
quantity  of  carbon  dioxide  for  various  -wave  bands  of  a  parallel  beam  of 
radiation.  His  formulas  for  the  equivalent  thicloiess,  c,  of  carbon  dioxide, 
and  for  the  monochromatic  absorptivity  of  carbon  dioxide,  in  a  column  of  air 
of  thickness  Z,  pressure  P  (in  atmospheres)  and  temperature  T  (degrees 


absolute)  are: 


c  '  8.7Z  f>i/T  , 


(5.4) 


R--  l~l/b  +  nCcffl)j .  D.5) 


In  equations  (5.4)  and  (5«5)»  c  is  the  length  of  absorbing  path  in  cr.is  of 
carbon  dioxide  at  II.T.P.  and  is  based  on  air  having  a  normal  carbon  dioxide 
content  (the  partial  pressure  of  carbon  dioxide  in  the  air  novr  being  .0052 
atmospheres).  The  constants  n  and  x  are  given  for  different  limits  of  the 
intense  absorption  band  centered  at  15  microns. 

Elsasser  |^5j  states  that,  in  the  case  of  square-root  absorption,  a 
slab  of  thicloiess  u  is  mathematically  equivalent  for  diffuse  radiation  to 
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a  linear  colunn  of  length  1«78  u  for  beaned  radiation.  Tho  condition  of 
square-root  absorption  is  satisfied  for  slabs  of  air  of  thioloiess  4  to  100 
meters  according  to  Panofsky  [lOj  •  ^°^   thinner  layers,  the  factor  1.78 
should  be  increased  slightly;  but  the  error  introduced  for  such  layers  will 
not  be  great  if  the  factor  1.78  is  used.  Inasnuch  as  diffuse  radiation 
occurs  in  the  atmosphere,  slab  thiclaiesses  must  be  increased  by  the  factor 
1,78;  \7hich  amounts  to  an  increase  in  the  equivalent  thicloiess,  c,  of 
equation  (5«4)  by  the  factor  1.78 •  Utilizing  this  corrected  equivalent 
thicloiess  in  the  absorption  fonnula  (3«5)  allovrs  this  formula  to  be  applied 
to  slab  thicknesses  of  the  atmosphere. 

In  this  study,  the  rate  of  cloange  of  temperature  of  layers  of  air  near 
the  ground  due  to  carbon  dioxide  as  \7ell  as  vra.ter  vapor  is  to  be  determined. 
Purthermore,  since  the  Atmospheric  Radiation  Chart  of  Elsasser  is  to  be  used 
for  the  rate  of  temperature-change  resulting  from  vra.ter  vapor,  consideration 
must  be  given  to  his  treatment  of  carbon  dioxide.  Elsasser  \^5j   considers 
that  the  absorption  by  carbon  dioxide  occurs,  for  the  most  part,  in  the  range 
752  cm"l  to  584  cm"^,  that  is,  bet\'/een  IJ.p  and  17.1  microns.  He  further 
considers  that  "very  thin  layers  radiate  as  black  bodies"  and,  therefore, 
neglects  the  radiative  effects  of  carbon  dioxide  in  comparison  to  those  of 
■water  vapor,  for  problems  concerning  net  flux.  The  chart,  however,  does 
give  a  separate  representation  of  the  black  body  flvix  contained  v.'ithin  the 
strong  band  intoirval  13 «p  -   17  •!  microns.  Callendar  j^  5 J  gives  the  values 
n  -  .35  aJid  X  ■  .55  as  sviitable  for  determining  the  absorptivity  by  (5) 
within  a  band  of  limits  15*0  -  17.0  microns.  It  is  assinned  that  the  values 
of  monochromatic  absorptivity.  A,  determined  from  (p. 4)  and  (p. 5)  are  appro- 
priate for  the  carbon  dioxide  limits  of  Elsasser,  13.3  -  17*1  microns. 


(8) 


p,\ 


ICnowing  the  ground  tenperature  and  the  mean  tenperature  of  the 
layers  of  atmosphere  under  consideration,  the  flux  through  the  top 
bo\uidary  of  each  layer  can  be  computed.  For  example,  for  the  first  slab 
above  the  ground,  the  equation  of  radiative  transfer,  soiaetimes  called 
Schv/arzchild's  equation,  gives 

Tftiere  D  is  the  flux  through  the  top  boundary  of  the  layer,  k-^   the  mono- 
chromatic absorptivity  of  the  layer,  E\q  the  black  body  flux  in  the  band 
13 •>  -  17 •!  microns  at  the  temperature,  T,  of  the  ground  (or  the  flux  into 
the  lower  boxmdary)  and  E^]_  is  the  black  body  f lu:c  at  the  moan  temperature 
f  of  the  layer.  Strictly  speaking  Schivarzchild^s  equation  (5«6)  is  true 
for  monochromatic  intensity,  or  beamed  flux,  of  radiation.  Actually 
Elsasser  £  5j  integrates  the  Schvrarzchild  equation  in  order  to  obtain  the 
flux  emerging  from  an  arbitrary  slab.  He  finds  that  this  integration  in- 
volves  the  computation  of  integrals  EcjfX)  ~     J  P        ^  ^y      where 

Ossec  y,  values  of  v/hich  are  given  in  |_5j  •  However,  he  also  shows  that 
these  computations  imply  the  increase  of  the  absorbing  path  by  the  factor  1.78 
discussed  above.  It  is  therefore  assumed  that  Sohvmrzchild's  equation  holds 
for  fluxes  E,  E^q,  and  Ej^  with  the  value  of  absorptivity  A  corresponding  to 
the  absorbing  path  c  replaced  by  1.78  c,  a  result  which  Elsasser  justifies. 
For  slabs  above  the  first,  equation  (5.6)  also  permits  the  computation 
of  upv/ard  flux  E  through  the  top  bomidary.  Hovrever,  for  such  computations 
E»  must  be  taken  to  represent  the  flux  into  the  bottom  of  the  slab,  which 
vreis  ths.t   out  of  the  top  of  the  slab  immediately  belc.T,  Ev,  would  still 
represent  the  black  body  flux,  in  the  band  1^.5  -  17 •!  microns,  corresponding 
to  the  mean  temperature  of  the  slab  under  consideration. 


(9) 


-o 


In  order  to  apply  the  same  procedure  to  the  problem  of  ooK^puting 
fluxes  for  the  dovravrard  radiation  of  carbon  dioxide  as  for  the  upward 
radiation,  it  is  necessary  to  know  the  fl;ix  entering  the  top  boundary  of 
the  highest  layer •  The  top  boundary  is  288  feet  for  this  stitdy.  Through- 
out the  period  used  in  this  study,  whenever  a  ground  inversion  exists,  this 
height  is  very  near  the  top  of  the  inversion.  The  slope  of  the  sounding 
above  the  groxmd  inversion  during  this  period  is  quito  similar  to  that  of 
the  standard  atmosphere.  This  remark  is,  on  the  whole,  true  also  of  the 
slope  of  the  sounding  above  the  288  feet  level  during  the  hours  in  vfhich 
no  ground  inversion  is  present.  Therefore,  the  atmosphere  from  10  kilo- 
meters is  divided  into  slabs  of  approximately  100  meters  thiolaiess  and 
having  the  appropriate  ten^erature  according  to  the  standard  atmospheric 
lapse  rate.  In  treating  equation  (5»6),  slab  1  is  to  be  considered  the 
topmost  slab,  slab  2  that  immediately  below,  etc.  It  will  be  shown  below 
that  the  amount  of  flux  from  above  10  kilometers  produces  negligible 
effect  at  the  288  foot  level.  Let  flux  E^  enter  the  top  of  slab  1.  The 
flvixes  leaving  the  bottoms  of  the  next  several  slabs  are  given  by  the 
follov/ing  table. 

Slab      Flux  through  the  lower  boundary 

1  1  E^  +  a, 

2  ml  Eq  ■♦-  ma  -».  b,  (5.7) 

3  nml  E  -♦•  nma  +  nb  -f  c, 

o 

In  (5.7 )»  1»  T^t   n etc.  correspond  to  the  factors.   (l-A-j_),  (l-Ap)) 

(1-A?),  etc*,  respectively,  of  equation  (3«^);  a,  b,  c  correspond  to 

the  values  A-j^  ■^Al*  -^  ■^A2*  "^^  ^^^*  *'••  ®'^°*  °-^  equation  (5#^)« 
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Since,  using  equations  (5 .4)  and  (3«5)»  the  average  value  of  1,  m, 
n,  •••••  etc.  is  of  the  order  0*5  for  a  layer  of  thickness  100  neters, 
and  Eq  is  of  the  order  of  10  cal  cm"2  (3  hr«)"^,  it  tiims  out  that  the 
contribution  of  the  initial  flux  E^  is  negligible  in  treating  equation 
(5.6)  for  the  dovmward  directed  flux.  "iTe  may  therefore  arbitrarily  set 
E  s  0  in  determining  the  flux  through  the  bottom  of  the  slabs  and  still 
obtain  the  correct  result  at  288  feet. 

Finally,  since  we  now  may  determine  the  amounts  of  upward  and  dovni^mrd 
directed  flux  at  every  level,  one  can  compute  net-upv/ard,  net-dovravvard  and 
consequently  net-outvfard  directed  fltix,  Fp  -  F, ,  for  every  consecutive  pair 
of  levels.  This  deterraines  the  radiative  heating,  or  cooling,  of  the  slab 
due  to  carbon  dioxide  radiation,  and  the  actual  rate  of  temperature-change 
can  be  computed,  as  before,  from  equation  (3«5)« 
5.   Heat  Transfer  by  lloleculax  Conductivity. 

The  upward  flux  of  heat  due  to  molecular  conduction  is  given  by  the 
formula: 

F«  -  -<<  dZ   ;  (5-8) 

where  oC  is  the  coefficient  of  molecular  conductivity  of  air,  and  has  the 
value  a  s  549(10)"^  cals  cm"-''  sec"*-^  T"-^  at  0°C,  but  is  a  function  of  tem- 
perature (c.f.  for  example,  p.  44  of  [ij  ). 

It  is  necessary  noxr  only  to  compute  upmrd  fluxes  F^^^  and  F^  for  the 
various  slabs  for  -which  lapse  rates  arc  knovjn,  and  again  use  equation  (5  •5) 
to  convert  net  flux  out  of  a  slab  into  a  corresponding  rate  of  change  of 
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temperature.  In  actual  practice  it  is  normally  found,  since  oc  is  very 
small,  tliat  in  order  to  register  a  measurable  net  flux,  the  changes  in 
lapse  rates  along  the  vertical  must  be  quite  extreme,  and  this  can  be 
expected  only  vrithin  the  lovrest  trro  or  three  slabs  • 
4«   Rate  of  Cliange  of  Temperature  by  Eddy  Conductivity. 

The  formula  for  the  observed  rate  of  temperature  change  is: 


(5.9) 


ar  wii)  ^  /Jit) + fix):  +f^ 

Zt        (*^/w      (.3t/c    U^</»n     \d^k) 
T/here  x'd'C  /yi^^^  due  to  T/ater  vapor, 

\d^  Jc  is  d^®  ^'^  carbon  dioxide, 
VB-T/^is  due  to  molecular  conduction, 

/2Z) 

\.dX  J^is   duo  to  eddy  conduction^ 


The  observed  rate  of  cliange  of  temperature  per  Ixalf-hour  is  given  for 
each  layer  by  the  data  used  in  this  study#  As  a  result  of  methods  previously 
described,  the  only  unknoi-m  of  fonmila  (p«0)  is  the  rate  of  change  of  tempera- 
ture due  to  eddy  conduction* 
5«   The  Coefficient  of  Eddy  Conduction. 

The  net-out\mrd  flux  for  a  layer  of  atmosphere  due  to  turbulent  con- 
duction may  be  expressed  as: 
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In  (5»10)  Kp  and  K-i  are  ohe  coefficients  of  eddy  conductivity  at  the  top 
and  bottom  of  any  slab  of  atmosphere,  respecti-vely.  Hiaking  the  substitu- 


tions ^  _:.  p^_ 


and 


^■^^¥9%^"^  > 


both  approximations  being  valid  to  vrithin  1%  for  the  layers  betvreen  the 
surface  and  288  feet,  we  obtain: 


In  (5»11),  all  items  labeled  with  subscript  1  refer  to  measurements  made 
at  the  center  of  slab  1,  all  items  with  subscript  2  refer  to  the  center  of 
slab  2j  j^  is  the  mean  of  p-|_  and  pp*  Thus  (5 .11)  determines  the  rate  of 
change  of  temperature  approximately  at  the  boundary  of  slab  1  and  slab  2« 
In  general,  slab  1  vdll  denote  any  lower  slab  and  slab  2  the  slab  immed- 
iately above  slab  1#  Thus  (5»ll)  nay  be  used  as  a  recursion  formula  to 
determine  all  values  of  Kj^  after  K^.i  has  been  determined.  The  only  limi- 
tation on  its  use  arises  when  -  F,^  =•  ^   ,  tlmt  is,  when  the  lapse  rate 
of  the  slab  is  dry  adiabatic ,  or,  not  significantly  different  from  Oj  • 

In  using  (5«ll)  as  a  recursion  formula,  it  is  necessary  to  determine 
a  value  of  "S.  at  some  specific  level  to  loermit  determination  of  all  other  Kjj_» 
In  order  to  do  this  without  making  too  many  assumptions,  a  value  of  K  at 
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sone  standard  leirel  can  be  determined  by  reference  to  the  theory  of 
nomentum-transfer  within  the  surface  la^'or*  The  coefficient  of  eddy 
viscosity  can  be  determined  at  a  standard  level,  say  16  feet,  and  it  can 
be  asstmed  to  give  a  representative  value  of  the  coefficient  of  eddy  con- 
diictivity.  The  particular  values  of  K  vfill  then  depend  to  some  extent  on 
momentum-transfer  theory,  but  the  comparative  values  reflect  only  the 
measurements  of  eddy  conductivity.  Recent  results  of  Lettatu  ["9 J  are 
drarm  upon  for  the  purpose  of  finding  a  value  of  E  at  sone  standard  level* 
According  to  the  Prandtl  theory  of  momentum-transfer*,  one  obtains: 


K-nm. 


(5.12) 


In  (5 •is)  1  is  the  moan  mixing  length  and  ^^   is  the  obsei^ed  shear  of 
the  horizontal  v/ind,  both  measured  at  a  given  level.  Lettau  f^l  gives  the 
folloTring  formulas  vdiich  permit  the  calculation  of  the  coefficient  of  eddy 
viscosity  Trithin  the  surface  layer: 


vrhere 


^I  "  "^  15T^   •  (5*14) 


♦HaurvTitz,  B.  Dynaraic  L^teorology.  Equation   (74.2).  HcGraw-Hill  1^41 . 
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In  (J.lp)  and  (p«14),  tho  sultscript  "a"  indicates  an  isotropic  or  adiabatio 
surface  la^^^r.   In  ^^  "^^^   mean  temperature  of  the  slab  for  vrhich.   ^^   '^^^ 

been  determined.  The  value  of  s-  may  be  detemined  by  the  Trell-loiovm 

logarithnic  lav/ 

sr  r  jz-i-z. 


«.-  T  f 


yo    '^"\    Z. 


(5.15) 


=  fl^    i,(ZfZ.) 


l^-y;;o-     KZ  +  ZJ     .  (3.16) 


In  (3«15)  and  (3»lo),  all  variables  are  understood  to  refer  to  the  surface 
layer*  For  example, 

t^  is  the  constant  value  of  the  stress, 

Ug_  is  the  horizontal  wind  speed, 

1  is  the  mixing  length, 

k  is  the  Yon  Ilarroan  constant,  k  -  .45, 

Z  is  the  roughness  parameter. 

PfT  ^ 

The  values  of  Zq  and  ■w — 9 —  can  be  obtained  from  a  plot  of   ^7    v/hen 
adiabatic  conditions  are  found.  Observations  of  Lettau  [8 J  indicate  that 
such  conditions  exist  at  sunrise  and  sunset.  An  average  Ug^  profile,  there- 
fore, can  be  dra-wn  based  on  both  sunrise  and  sunset  data  in  order  to  deter- 


mine Z  and  the  mean  value  of  "Vj  — ^^ —    for  the  day. 


Ik 
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IV,  RM:ES  of  CZAilGS  OF  TEIPERATURE 

Im      Rate  of  Change  of  Temperature  Due  to  Tfater  Vapor. 

Radiation  transfer  due  to  vraiter  vapor  T/as  computed  for  each  la^ner  of 
the  soundings  chosen  from  the  data  "by  the  method  described  in  Chapter  III, 
Section  1,  of  this  study.  This  method  requires  a  laaov/ledgo  of  the  specific 
humidity  at  each  level.  The  data  [6j  gives  the  vapor  pressure  at  each 
level.  These  vapor  pressures  -were  converted  to  specific  Iiunidity  by  use 
of  the  Kiefer  I5alti-Prossure  Hyr;rometric  Chart.  Since  equation  (5.1)  re- 
quires the  pressure-thiclaiess,  ^p,  of  all  layers  it  Tms  necessary  to  com- 
pute the  pressure  at  every  level  by  equation  (3.2).  To  do  this,  the  sur- 
face pressure  must  be  knoTm.  The  surface  pressure  vras  computed  from  the 
sea  level  pressure  recorded  for  Austin,  Texas  on  the  surface  map  at  the 
middle  of  the  period  under  study,  27-28  September  1948.  This  Tra.s  done  by 
reversing  the  procedure  of  the  reduction  of  pressures  to  sea  level,  knovring 
the  height  of  the  base  of  the  totTer  (525  feet  IJSL). 

A  (u,T)  relationship  y/as  computed  for  all  levels  of  the  soundings  and 
plotted  on  the  Slsasser  Cliart.  The  difference  of  the  net  fluxes  of  each 
slab  of  the  soundings  Y/as  obtained  by  making  the  required  area  measurements 
with  the  aid  of  a  planimetcr.  Figures  1  and  2  show  schematically  the  areas 
that  represent  the  differences  of  net  flux,  for  one  slab  of  atmosphere  for 
the  folloT/ing  cases: 

(a)  The  slab  lies  v;ithin  a  ground  inversion  (Figure  l). 

(b)  The  temperature  decreases  vri.th  height  (Figure  2). 

In  Figures  1  and  2,  the  reference  levels  L,  and  Lp  denote  respectively  the 
bottom  and  top  of  the  slab,  both  tal:en  to  lie  belovr  the  288  foot  level. 
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Difference  of  llet  Flvixos  for  a  Slab  Belov/  288  Feet 

Figure  1 
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Consider  next  Figure  1.  The  (u,T)  oxarre   extending  above  L-j_  must 
evon-tiially  cross  that  above  L2  because  the  fomer  curve  eventually  corres- 
ponds to  a  greater  integrated  optical  depth  of  water  vapor.  Suppose  this 
"cross-over"  point  is  labeled  A,  At  2200  LCT,  the  complete  radiosonde  T,ms 
available  to  supplement  the  2200  nicrometeorological  sounding  •  It  -.-as  found 
that  the  cross-over  points  for  the  (u,x)  plots  above  the  slab  levels  Lj   and 
Lp  occurred  from  about  5OO  feet  for  the  loTrer,  thinner  slabs  to  about  8 00 
feet  for  the  higher,  thicker  slabs.  Ijoreover,  it  vreis  found  that  the  area 
computation  made  by  terminating  the  (u,T)  curves  above  L^  and  L2  at  points 
Bt  and  B2  corresponding  to  the  288  feet  level  Tra.s  in  error  only  by  2-1% 
•when  coii^ared  vd-th  the  areas  given  by  the  complete  (u,T)  curves*  This  small 
error  is  due  to  a  tendency  of  positive  and  negative  areas  beloiT  the  isotherm 
B^Bo  to  cancel  each  other* 

Consider  next  Figure  2,  representing  typical  daytime  conditions.  IJo 
exact  comparison  is  laiowi  here  of  the  error  made  in  neglecting  the  area 
below  and  to  the  right  of  the  isotherm  B]^B2«  Ho  cross-over  point  occurs  in 
this  case  so  there  is  no  tendency  for  cancellation  of  positive  and  negative 
areas.  For  the  lower,  thinner  slabs  little  error  could  occur  because  the 
(u,T)  curves  above  the  slab  boundaries  are  in  general  nearly  coincident  at 
B-i  and  B^*  The  greatest  error  should  occur  in  the  slab  185-255  feet  and  the 
error  should  be  progressively  sxaaller  in  the  slabs  immediately  belovf.  The 
neglected  area,  if  appreciable,  should  cause  a  reduced  amount  of  cooling  or 
increased  vrarming  in  the  midday  hours  in  the  slab  185-255  feet.  Actually 
Table  1  shoT;s  no  evidence  of  reduced  cooling  of  the  slab  185-255  feet.  For 
tho  ratio  of  the  I3OO  LCT  cooling  of  this  slab  to  those  of  the  threo  slabs 
immediately  below  is  almost  identical  with  the  corresponding  slab  ratios  of 
the  2200  LCT  sounding,  whose  values  are  baown  to  be  correct. 


(19) 


I    0 
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(Times  given  in  LOT.) 


Layer 
(in  feet) 

0200 

AF 

0500 
AF 

0800  ' 

AF 

1100 
AF 

1300 

AF 

1600 
AF 

1900 

AF 

2200 

AF 

0-,5 

.01 

-.01 

.03 

.04 

•02 

.00 

-.045 

.00 

.5-1.5 

.045 

•05 

-.04 

-.40 

-.38 

-.175 

-.02 

.04 

1.5-3.0 

.065 

.065 

-.10 

-.44 

-.48 

-.21 

.03 

.04 

3-6 

.08 

.08 

-.22 

-.56 

-.51 

-.28 

.05 

.10 

6-12 

.10 

•10 

-.23 

-.87 

-.67 

-38 

.11 

.08 

12-20 

.08 

.11 

-.17 

-.66 

-.75 

-.38 

.19 

.08 

20-35 

.10 

.10 

-.21 

-.80 

-.79 

-.34 

.23 

•08 

35-55 

.08 

•10 

-.24 

-.82 

-.79 

-.35 

•33 

•11 

55-80 

•02 

.08 

-.29 

-.66 

-.56 

-.32 

.26 

.19 

80-110 

.04 

.07 

-.02 

-.61 

-•60 

-•34 

•28 

•52 

110-145 

.13 

.05 

-.07 

-.54 

-.59 

-•31 

.21 

.36 

145-185 

•19 

•25 

-.33 

-.56 

-.57 

-•28 

.20 

•29 

185-235 

.21 

.26 

-.17 

-•45 

-.46 

-.28 

.20 

•28 

Heat  Transfer  Due  to  Tfater  Vapor,     cal^   (3  hr.)"^ 


TABLE  1. 
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(Times  given  in  LCT.) 


Layer 
(in  feot) 

0200 

0500 
(°c) 

o8oo 

(OC) 

1100 
(Oc) 

1500 

(Oc) 

1600 
(OC) 

1900 
(oc) 

2200 

(°9) 

0-.5 

-2.28 

2.28 

-b.op 

-  9. $5 

-  5.47 

.00 

10.25 

.00 

.5-l«5 

-5.12 

-5.54 

4.56 

46.86 

59.85 

20.50 

2.28 

-4.56 

1.5-5.0 

-4.94 

-4.84 

6.41 

54.04 

57.15 

16.24 

-2.28 

-5.04 

3.0-6 

-5.01 

-5.01 

8.55 

21.87 

19.91 

10.95 

-1.92 

-5.80 

6-12 

-1.88 

-1.88 

4.59 

16.90 

15.08 

7.55 

-2.11 

-1.52 

12-20 

-1.15 

-1.55 

2.45 

9.65 

10.90 

5.56 

-2.74 

-1.14 

20-55 

-  .75 

-  .75 

1.59 

6.18 

6,20 

2.64 

-1.76 

-  .61 

55-55 

-  .44 

-  .56 

1.57 

4.77 

4.59 

2.04 

-1.91 

-  .61 

55-80 

-  .09 

-  .56 

1.52 

5.08 

2.61 

1.50 

-1.21 

-  .88 

80-110 

-  .15 

-  .26 

.76 

2.55 

2.58 

1.51 

-1.08 

-1.24 

110-145 

-  .45 

-  .29 

.25 

1.8o 

1.97 

1.04 

-  .70 

-1.20 

145-185 

-  .54 

-  .70 

•95 

1.75 

1.65 

.81 

-  .58 

-  .85 

185-255 

-  .52 

-  .65 

.41 

1.05 

1.14 

.70 

-  .50 

-  .69 

Three-hourly  Hate  of  Change  of  Temper atiire  Due  to  'Jater  Vapor 
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Table  1  gives  the  values  obtained  for  heat  transfer  due  to  vrciter 
vapor  in  cal.  (5  hr.)"  •  Table  2  gives  the  three-hourly  rate  of  change 
of  temperature  for  each  layer  of  each  soxinding,  computed  from  Table  1  by 
use  of  equation  (5 •J)* 
2*   Rate  of  Change  of  Temperature  Due  to  Carbon  Dioxide. 

The  differences  in  upvjard-direoted  flux  due  to  the  radiative  transfer 
of  carbon  dioxide,  in  the  band  15 •?  -  17«1  microns ,v;ere  computed  for  each 
layer  by  the  use  of  equation  (5»b),  The  differences  in  dovnav/ard-directed 
flux -were obtained  from  the  method  of  Chapter  III,  Section  2« 

The  computations  of  dovravvard-directed  flux  took  into  account  that  the 
mid-period  sounding  vreis,  on  an  average,  ten  degrees  vfarmer  throughout  than 
that  of  the  standard  atmosphere,  although  the  slope  vjas  assimed  to  be  that 
of  the  standard  atmosphere.  This  amounted  to  increasing  the  average  emission 
from  the  layers  by  the  factor  1(^,  ■which  factor  is  the  percentage  increase 
in  flux  Evj^  of  equation  (5«^)  due  to  the  increase  in  the  moan  temperature, 
Tjj^,  of  each  layer  compared  to  the  standard  atmospheric  temperature.  The 
adjusted  value  obtained  for  the  flux  through  the  top  boimdary  of  the  layer 
288-255  feet  Tras  18.79  cal  cm"^  (5  hr.)"^.  This  value  -was  used  as  a  constant 
value,  for  the  period  of  this  study  and  computations  of  flux  differences  for 
the  layers  of  each  soxmding  under  study  v/ere  made  using  equation  (5.6).  The 
values  of  black  body  flux  of  emission  for  the  carbon  dioxide  band  15 "5  -  17«1 
microns,  at  the  mean  temperature  of  oach  slab  vrere  taken  from  the  top  of  the 
Atmospheric  Radiation  Cliart.  The  values  of  monochromatic  absorptivity  -ivere 
computed  using  equation  (5«4)  and  (5»5)» 
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9"  jo'xotj: 


0-:^ 


'L^-L^i 


..A 


(Times  given  in  LOT.) 


Layer 
(in  feet) 

0200 
AF 

0500 
AF 

0800 

AF 

1100 
^F 

1300 
AF 

1600 

i^F 

1900 

2200 

0-.5 

-•012 

-.010 

-.00 

-.010 

-.010 

-.003 

-.025 

-.010 

.5-1.5 

-•005 

-.007 

-.015 

-.064 

-.062 

-.030 

-.021 

-.005 

1.5-3.0 

-.001 

-.002 

-.027 

-.101 

-.094 

-.044 

-.007 

-.005 

3.0-6.0 

.004 

.003 

-.042 

-.149 

-•133 

-.062 

.010 

.002 

6-12 

.015 

.016 

-.045- 

-.211 

-.187 

-.083 

.048 

.015 

12-20 

.012 

.013 

-.053 

-.216 

-.222 

-.089 

.074 

.018 

20-35 

.012 

.008 

-.080 

-.269 

-.285 

-.113 

.113 

.031 

35-55 

.022 

.010 

-.095 

-.368 

-.286 

-.114 

.155 

.064 

55-80 

.023 

.005 

-.102 

-.232 

-.257 

-.106 

.169 

.122 

80-110 

.044 

.006 

-.066 

-.207 

-.241 

-.092 

.174 

.185 

110-145 

.092 

•023 

-.082 

-.184 

-.224 

-.056 

.184 

.208 

145-185 

.127 

.051 

-.124 

-.115 

-.159 

-.003 

.210 

.212 

185-235 

.165 

.092 

-.126 

-.037 

-.087 

.061 

.24G 

.231 

Heat  Transfer  Due   to  Carbon  Dioxide,  .Cal.(3  hr.) 
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(Times  ^iven  in  LOT.) 


Layer 
(in  feet) 

0200 

(°c) 

0500 
(°c) 

o8oo 

(OC) 

1100 

(°c) 

13C0 

(OC) 

l6oo 

(OC) 

1900 

(OC) 

2200 
(Oc) 

0-.5 

2.71 

2.30 

.16 

2.41 

2.84 

.77 

5.60 

2.26 

.5-1.5 

•62 

.73 

1.66 

7.52 

7.30 

3.47 

2.42 

.55 

1.5-5.0 

.11 

.17 

2.07 

7.81 

7.25 

5.45 

.52 

.37 

5.0-6 

-.14 

-.12 

1.57 

5.84 

5.20 

2.41 

-.57 

-.06 

6-12 

-.28 

-.30 

.93 

4.10 

3.65 

i.6o 

-.93 

-.29 

12-20 

-.17 

-.18 

.76 

3.15 

3.23 

1.51 

-1.07 

-.25 

20-35 

-.08 

-.06 

.60 

2.08 

2.24 

.87 

-.87 

-.24 

35-55 

-.12 

-.06 

.54 

2.13 

1.66 

.67 

-.90 

-.35 

55-60 

-•10 

-.03 

.47 

1.08 

1.19 

.50 

-.78 

-.56 

80-110 

-.16 

-.03 

.25 

.80 

.95 

.55 

-.66 

-.71 

110-145 

-.30 

-.08 

.27 

.62 

.75 

.19 

-.62 

-.69 

145-165 

-.36 

-.14 

•35 

.36 

.45 

.00 

-.61 

-.61 

185-235 

-.35 

-.23 

.30 

.08 

.22 

-.15 

-.62 

-.57 

Three-hourly  Rate  of  Temper  attire-change  Due  to  Carbon  Dior^ide 
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Table  3  gives  the  results  obtained  for  net-out\mrcl  radiation* 
Tabic  4  rives  the  resultant  rate  of  temperattare-change  for  each  slab 
of  each  sounding,  computed  using  the  results  of  Table  3  and  equation  (3«3)« 
3«   Resxiltant  Radiative  Rate  of  Temperature-Change* 

The  tenperature-change  values  of  Tables  2  and  4  -were  added  algebraically 
yielding  a  resultant  radiative  teraperaturo-change  for  the  slabs*  Since 
these  tenperature-clionges  are  to  bo  compared  with  the  observed  teinperature- 
chjange  at  the  various  boundeiries  it  -was  necessary  to  adjust  them  to  these 
hoimdaries*  This  ■was  done  by  considering  the  radiative  temperature-changes 
as  existing  at  the  midpoints  of  the  slabs  and  then  ijiterpolating  to  the 
slab  boxmdaries  on  log-log  paper*  These  results  are  presented  in  Table  5* 
wherein  all  three-hourly  rates  have  been  reduced  to  half-hourly  rates  by 
dividing  by  6* 
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(Times  given  in  LCT) 


Bovtadary 
(in  feet) 

0200 
(OC) 

0500 
(°c) 

o8oo 

(OC) 

1100 
(OC) 

1300 
(OC) 

1600 

1900 
(OC) 

2200 
(OC) 

.5 

-•20 

.24 

-.40 

2.21 

2.77 

1.42 

2.02 

.03 

1.5 

-.77 

-.79 

1.19 

8.23 

7.68 

3.71 

.35 

-.58 

3.0 

-.71 

-.69 

1.49 

6.19 

6.34 

2.95 

-.32 

-.51 

6.0 

-.47 

-.47 

1.40 

4.25 

3.72 

1.98 

-.42 

-.53 

12 

-.30 

-.33 

.73 

2.88 

2.60 

1.33 

-.57 

-.27 

20 

-.19 

-.24 

.47 

1.87 

2.06 

.95 

-.57 

-.20 

35 

-.12 

-.12 

.34 

1.28 

1.25 

.52 

-.45 

-.15 

55 

-.06 

-.09 

.31 

.?5 

.86 

.40 

-.41 

-.19 

80 

-.04 

-.06 

.23 

.60 

.59 

.30 

-.31 

-.28 

110 

-.08 

-.05 

.13 

.47 

.50 

.24 

-.26 

-.32 

145 

-.13 

-.10 

.14 

.38 

.41 

.17 

-.21 

-.28 

185 

-.15 

-.14 

.17 

.28 

.30 

.12 

-.20 

-.23 

235 

-.15 

-.11 

.10 

.16 

.20 

.06 

-.18 

-.19 

Half-hotirly  Rate  of  Temperature -change:  V/ater  Vapor  +  Carbon  Dioxide 


TABIJE  5. 


(26) 


•n*- 


4«   Rate  of  Change  of  Temperature  Dv.e  to  Molecular  Conduotion. 

The  flux  due  to  nolectxlar  conduction  vrcis  computed  by  the  use  of 
equation  (5»8)«  The  temperature,  riven  in  the  observed  data,  for  each 
level  -was  used,  hence  the  fliix  -was  considered  as  passing;  through  the 
center  of  each  layer.  Consequently,  the  temperature-change  computed 
from  the  value  of  the  flux  differences  for  tvro  consecutive  layers  must 
be  ccisidorod  as  beinr  at  the  boundary  betvreen  the  two  layers.  This 
resulted  in  a  rate  of  temperature-change  at  the  boundary  of  each  slab 
vd-thout  the  necessity  of  interpolation.  The  results  obtained  for  rate 
of  temperatixre-change  due  to  molecular  conduction  are  given  in  Table  6. 

(Times  given  in  LCT) 


Boundary 
(in  feet) 

0200 
(OC) 

0500 
(OC) 

o8oo 

(OC) 

1100 
(OC) 

1500 
(°c) 

1600 

(OC) 

1900 
(OC) 

2200 

(OC) 

.5 

-.90 

-1.90 

-i.6o 

7.20 

8.0 

4.10 

-.57 

-2.60 

1.5 

-.12 

-.16 

-.09 

*51 

.45 

.17 

-.27 

0.0 

3.0 

.06 

.07 

.05 

-.06 

Half -hourly  Rate  of  Temperature-change  Due  to 
Molecular  Conduction 


TABLE  6. 

5.   Rate  of  Change  of  Temperature  Due  to  Eddy  Conduction. 

These  results  follow  innnediately  by  the  use  of  the  observed  ^;^ 
and  the  values  of  Table  5  and  6  in  equation  (3.9).  The  values  for  the 
rate  of  temperature-change  due  to  eddy  conduction  are  given  in  Table  7* 


ar 


(27) 


(Tines  gi-ven  in  LCT.) 


Boundary 
(in  feet) 

0200 
(oc) 

0500 

0800 

1100 
(OC) 

1500 

(°c) 

1600 
(OC) 

1900 

(^c) 

2200 
(OC) 

•5 

.80 

1.66 

1.60 

-9.01 

-10.57 

-6,52 

-1.05 

-1.67 

1.5 

0*59 

0.75 

1,22 

-7.60 

-7.65 

-4.48 

-.02 

0.18 

5 

0.21 

0.49 

0.47 

-4.95 

-5.91 

-5.48 

-.02 

0.01 

6 

-0.25 

0.17 

0.50 

-5.15 

-5.62 

-2.28 

-0.28 

-0.07 

12 

-0,50 

0.05 

0.97 

-2.18 

-2.20 

-1.25 

-0.55 

-0.25 

20 

-0.51 

0.04 

1.45 

-1.57 

-2.56 

-1.15 

-0.05 

-0.50 

35 

-0.48 

-0.28 

1.46 

-0.78 

-0.85 

-0.52 

0.15 

-0.55 

55 

-0.54 

-0.11 

1.49 

-0,45 

-0.56 

-0.50 

0.51 

-0.21 

80 

-0.06 

-0.04 

1.57 

-0.00 

-0.59 

-0.10 

0.21 

-0.22 

110 

0.28 

0.05 

0.47 

-0.07 

-0.10 

-0.14 

0.06 

-0.28 

145 

0.05 

0.50 

1.46 

0.02 

•  -0.01 

0.05 

0.01 

0.28 

185 

-0.05 

0.14 

1.55 

0.02 

0.20 

0.08 

-0.10 

0.45 

2^5 

-0.55 

0.12 

1.50 

•04 

0.10 

0.04 

0.08 

.59 

Half-hourly  Temperatiire-change  Due  to 
Eddy  Conduction 
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V,  RESULTS  AIID  COIICLUSIOIIS 

1.   TTothod  of  Computation  of  the  Coefficient  of  Eddy  Conductivity. 

As  indicated  in  Chapter  III,  Section  5*  it  "vms  decided  to  determine 
for  all  soundings  evaluated  a  first  value  of  K  from  the  theory  of  turbu- 
lent transfer  of  momentina*  For  this  rMiT'nnRA  it  vra.s  necessary  to  maJ-re  use 
of  observed  -vrinds  in  order  to  evaluate  ^Z        ^o^  ^s®  ^   equation  (3 -IS). 
TTfo  vrinds  yrere   riven  consistently  in  the  data  vfhich  vrere  alivays  vrithin  the 
surface  layer,  one  at  12  feet  and  one  at  41  feet.  For  each  of  the  ei^ht 
soundings  studied,  these  tvrt)  values  of  the  wind  were  plotted  against  log  Z 

on  semi-log  paper,  and  a  straight  line  joining  the  tvro  points  was  drawn. 

dg 
The  value  of  ^Z   ^Ji  the  vicinity  of  l6  feet  was  then  read  off  as  being 

Au 

the  value  /^Z    corresponding  to  a  height  increment  of  100  centimeters 

at  l6  feet.  The  reasons  for  choosing  the  l6  foot  level  as  the  standard 
reference  level  are  tiTO-fold: 

(a)  The  level  had  to  be  within  the  lov/est  levels  for  which  wind 
data  are  available. 

(b)  The  level  had  to  be  within  the  lowest  slab  consistent  vfith  (a) 
for  v/hioh  lapse  rate  data  is  available.  This  is  necessary  since,  in 
equation  (5»15)»  "tbe  mixing  length,  1,  involves  di'Z    *  V.dZ     Qy   • 

In  satisfying  requirement  (b),  it  ims  necessary  to  use  the  slab  12-20 

de 

feet,  Y/hich  vreis  considered  to  define  dZ    at  l6  feet. 

In  several  cases  however,  it  was  necessary/  to  shift  to  the  level  23«5 
feet.  This  necessity  arose  when  an  obviously  inconsistent  lapse  rate  was 
indicated  in  the  layer  12-20  feet.  This  shifting  to  the  level  2^.5  feet 
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amounted  to  dealing  with  tlie  over-all  lapse  rate  in  the  layer  12-55  feet, 
v/hich  then  gave  ^  AT    at  25,5  feet. 

In  order  to  evaluate  the  nixing  length,  1  (the  remaining  i«rn  of 
equation  (5«12))  at  the  standard  level,  a  knowledge  of  the  corresponding 
value  Ij^,  the  mixing  length  in  an  isotropic  atmosphere,  is  necessary^ 
According  to  Lettau  [9]  isotropic  or  adiabatic  turbulence  is  most  nearly- 
realized  arovmd  sunset  or  sunrise*  Thus,  in  order  to  derive  the  adiabatic 
velocity-profile  for  the  24-hour  period  under  investigation,  an  average  of 
three  half-hourly  Trind-speed  reports  at  both  sunrise  and  sunset  was  obtained. 
Then,  plotting  mean  wind-speed  against  log  Z  on  semi-log  paper  for  the  tvro 
values  of  Z  a  straight  line  vra.s  obtained.  Figure  p,   which  according  to 
Lettau  should  represent  the  wiiid  profile  for  Iianor,  Texas  under  adiabatic 
conditions  for  27-28  September  1948*  The  straight  line  obtained  apparently 
does  represent  a  close  approximation  to  the  true  adiabatic  atmosphere,  since 
Figure  J>   indicates  a  roughness  parameter  Z  -  2  cms,  which  checks  rather  vrell 
with  values  quoted  by  Stewart  [  1 J  for  fallow  fields.  Knovring  the  value  of 
the  roughness  parameter  permits  the  computation  of  1  from  equation  (5»13)» 

The  value  of  a^    at  16  foot  -vvas  then  read  directly  off  the  velocity 
profile  of  Figure  5»  Since  I6  feet  4  4^0  cms  this  can  be  done  quite  simply 
by  reading  off  the  wind-speed  differences  bet^veen  60O  and  400  cms  and  divid- 
ing by  200  cms*  This  calculation  leads  to  •15  sec"'^,  a  value  v;hich  may  be 
checked  mathematically  using  the  theoretical  equation  of  the  straight  line 
of  Figure  5»  This  is  done  below. 
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Figure  5« 
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The  equation  of  the  wind  profile  line  of  Figure  5  ccin  be  v<Titten 
as  in  (5,15)«Using  the  nean  value  of  Ug_  at  the  12  foot  level  (566  cms), 
k  w  .^ 
Then  from  equation  (I6), 


►45  J  Z  ■  2  cms,  Ug^  «  5«8  mps,  vre   obtain  "\ I -^-  -  52 #8  cm  sec"^. 


(■ 


r^l      =-:4r"^52.8)  i-^)  z   .15  seo-1. 

Similarly  at  25«5  feet  the  value  of  ^V         ""^.s  fotmd  to  be: 

^     =^(52.8)^=.10sec-l. 

At  these  ti-ro  levels  Ig^  s  220  and  523  cms,  respectively.  Values  of  K  at 
the  reference  level  may  now  be  determined  once  the  values  of   ^  Z  ■>  ^c 
^2  "^Q  Imotm.  These  results  are  tabulated  below  in  Table  8  for  the 
various  tines  of  day: 


02C0* 

0500* 

1100* 

1500 

IcOO 

1900 

2200 

de 

3Z 

"y^*^   -OOOo? 

.00058 

-.00075 

-.0028 

-.00061 

.00092 

•00092 

e. 

'^^cm  .0050 

•0050 

.0031 

.0066 

.0030 

.0066 

•0065 

sec''  .185 

.54 

•215 

.458 

•222 

•214 

•438 

K     c-mse^^^*^^^         28,000     40,900     63,800    38,000    8,000     16,000 


Computation  Values  for  Determination  of  K  at  Reference  Level 


TABLE  8^ 
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In  Tablo  8,  the  *  alongside  the  time  indicates  tha..  >.:.o-  rGf'^->''^--""3 
level  %ms  taken  at  23*5  feet  due  to  the  non-representative  lapse  rates 
betT.'reen  12  and  20  feet*  The  0800  values  of  K  were  not  computed  since 
many  of  the  slabs  had  either  non-representative  or  dry  adiabatic  lapse 
rates  so  that  the  heat  transport  equations  (5»11)  yield  no  solutions. 

For  all  soundinrs,  other  than  the  0800  sounding,  it  was  possible  to 
obtain  solutions  of  equations  (3«ll)  using  the  values  of  K  of  Table  8  as 
a  starting  point.  The  values  of  K  for  0200,  0500,  1100,  1300,  l600,  1^00 
and  2200  LCT  at  the  indicated  elevations  are  listed  in  Tablo  9»  s-nd  shorm 
graphically  in  Plate  I»  In  Table  9  the  occurrence  of  blank  spaces  indi- 
cates the  existence  of  a  non-representative  lapse  rate  in  the  slab  having 
the  indicated  height  as  its  midpoint.  It  vra,s  therefore  necessary  to  con- 
solidate adjacent  slabs,  obtaining  a  value  of  K  corresponding  to  the  new 
midpoint,  as  discussed  earlier  in  Section  1  of  this  chapter. 


♦Actually,  isothermal  lapse  rates  vrere  not  significantly  different  from 
the  dry  adiabatic  for  many  of  the  laj/ers  under  study,  when  the  gust  error 
of  "t   #2^0  is  considsrod* 
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(Times  given  in  LOT.) 

Level     0200      0500      1100      IJOO      loOO      1^00      2200 
(feet)   I[(10)^    K(10)_^_   hC^P)?.  .  .  ^(lO)^    K(lO)?    K(10)^    K(10)^ 


.25 

.075 

.14 

.045 

.26 

.04 

.056 

.09: 

1.0 

.51 

.68 

.50 

5.7 

.06 

.16 

2.25 

1.9 

5.0 

2.4 

16.4 

1.7 

.41 

2.5 

4.5 

1.5 

2.4 

5.8 

180. 

.75 

2.5 

7.5 

15. 

9.0 

4.9 

5.9 

20. 

48. 

1.5 

8.7 

16,0 

64. 

58. 

8.0 

16. 

25.5 

15. 

28. 

41. 

27.5 

500. 

65. 

5.1 

9.2 

45.0 

11. 

25. 

460. 

1700. 

55. 

8.0 

11. 

67.5 

17. 

28. 

500. 

150. 

52. 

9.0 

82.5 

120. 

95.0 

8. 

21. 

75. 

20. 

127.5 

10. 

18. 

550. 

25. 

550. 

165.0 

17. 

20. 

28. 

840. 

210.0 

17. 

82. 

210. 

920. 

261.5 

18. 

87. 

190. 

1500. 

The  Coefficient  of  Eddy  Conduction,  K,  (cn^  sao"*-^) 


TABLE  9, 


(M) 


In  order  to  show  moro  clearly  the  variations  of  K  with  both  height 
and  time,  graphs  of  K  at  fixed  elevations  were  draTm  showinr  the  tine- 
variation  of  K  at  these  elevations.  The  fixed  elevations  used  are  J 
inches,  6,  50^  ^0  and  125  feet  and  the  results  are  shoim  in  Figures  4, 
5,  6,  7,  and  8. 
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Diurnal  Variation  of  "K"  at  .25  Feet 

Figure  4. 
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Figure  5« 
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Figure  7 
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Figure  8# 
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2«   Discussion  of  Results, 

Lettau  [^yj  ^.las  define'  '"   surface  la^er  essentially  as  the  layer 
in  vfhich  the  turbulent  time-rate  of  change  of  a  property  is  noslirible 
in  comparison  v;ith  its  turbulent  flux.  In  the  solution  of  equations  (5*ll) 
for  the  various  levels,  it  is  evident  that  the  value  of  the  lapse  rates  be- 
come exceedingly  critical  as  the  value  of  the  lapse  rate  approaches  the  dry 
adiabatic  0^   •  For,  if  the  maximum  ;nAst  error  '^•2  ^C   occurred  in  opposite 
senses  at  adjacent  levels,  it  would  cause  an  error  in  observation  of  the 
lapse  rate  by  several  times  0,   ,  depending  upon  the  layer  thiolaiess.  Thus 
it  follows  that  the  nest  representative  values  of  K  may  be  expected  vrfiere 
the  lapse  rate  y  is  largest  compared  to  Jj  and  especially  vjhere  it  is  of 
the  order  of  50-100  (or  more)  times  dj  ♦  In  order  to  emphasize  this  de- 
pendence upon  lapse  rate,  let  us  consider  the  values  of  the  lapse  rate  in 
the  various  layers  for  the  0200  and  1^00  LOT  sottndin2;s»  These  are  shovm 
in  Table  10#  In  this  table,  the  abbreviation  "n»r«"  indicates  a  non- 
representative  lapse  rate;  therefore,  necessitating  a  consolidation  of 
the  two  adjacent  layers. 

Notice,  from  Table  10,  that  at  1^00  LCT  the  lapse  rates  were  strongly 
super-adiabatio  near  the  ground  but  not  significantly  different  from  the 
dry  adiabatic  above  80  feet#  At  0200  LCT  the  lapse  rates  vrere  more  nearly 
uniform  than  at  I5OO  LCT#  A  ooiq^arison  shows  that  the  IJOO  LCT  lapse  rates 
were  numerically  greater  than  those  for  0200  LCT  betx'reen  the  surface  and  6 
feet,  but,  above  6  feet,  the  early  morning  lapse  rates  in  general  exceed 
those  of  mid-day#  For  purposes  of  comparison,  the  value  of  the  dry  adiabatic. 
lapse  rate  is  Oj  z  •^o  x  10"^  °C  per  cm» 
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Lapse  Rate  (°c/cn)  x   10"4 


Laj'er 
(feet) 

0200  LCT 

1300  LCT 

0-.5 

-1 ,1x105 

7 .1x10^ 

.5-1.5 

-1 ,6x102 

4,9x10^ 

1.5-3 

-44 

110 

3-6 

-55 

11 

6-12 

-16 

38 

12-20 

n,r. 

29 

12-35 

-5.7 

13 

20-35 

n,r. 

4.4 

35-55 

-6,6 

4,9 

55-80 

-3.9 

1.3 

80-110 

-9.8 

8.0 

110-145 

-7.5 

isothermal 

145-185 

-4,1 

ti 

185-235 

-3.9 

R 

Values  of  Lapse  Rates  at  0200  and  I3OO  LCT 


TABLE  10, 


(42) 


In  general,  it  ■vjb.s   found  that  v/hen  Y  was  significantly  different 
from  3j  ,  the  turbulent  rate  of  change  term  of  equations  (5«ll)  "was 
conpletely  negligible*  Due  to  the  inability  to  smooth  out  the  turbulent 
fluctuations  of  teir^erature ,  the  definition  of  the  surface  layer  in  this     j 
study  is  the  layer  in  -which  the  lapse  rate  is  significantly  different  from 
dry  adiabatio,  both  the  gust  and  instrumental  errors  being  considered.  Thus 
the  surface  laj^er  at  I5OO  LCT  extends  to  about  80  feet,  but  at  0200  LOT  ex- 
tends to  at  least  235  feet.  This  is  evident  from  the  values  of  Table  10. 
It  must  be  realized  that  a  negative  value  of  this  lapse  rate  is  more  signifi- 
cantly different  from  Oj     than  a  corresponding  positive  value  of  the  lapse 
rate.  Then,  for  all  practical  purposes,  equation  (5.11)  vfithin  the  surface 
layer  reduces  to: 


k,  (If ' -i),  -  ff .  (-if  *  ^)r  ° . 


(5.1) 


■where  subscripts  1  and  2  refer  to  a  lower  and  upper  level,  respectively. 

Let  us  consider  the  actual  values  of  K  shovm  in  Plate  I  in  vrhich  the 
value  of  K  is  plotted  on  a* log-log  scale  against  the  height  of  the  midpoint 
of  the  slab  to  -wiiich  the  value  of  K  corresponds.  It  is  noted  that,  for  any 
sounding,  the  first  three  points  in  order  of  ascending  height  are  very  close 
to  a  straight  line,  there  being  some  deviation  from  the  straight  line  re- 
lationship as  the  fourth,  fifth,  etc.,  points  are  considered.  This  deviation, 
commencing  vrith  the  fourth  point,  is  due  mainly  to  two  factors: 

(1)  Decreasing,  therefore  less  representative,  values  oi  the  lapse  rate. 

(2)  The  change  in  temperature  elements  above  the  three  foot  level;  there 
being  non-aerated  resistors  below  and  aerated  resistors  above. 
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Therefore,  it  -vras  considered  pemissible  to  fit  the  "best  line  to  each 
sounding  in  Plate  I  up  to  the  naxiraum  height  for  vfhich  K  and  Z  both  in- 
crease together*  The  line  of  best  fit  for  each  sounding  vras  dravm  in 
accordance  vath  the  follovfing  principles: 

(a)  Through  the  first  point, 

(b)  Through  as  mainy  additional  points  as  possible. 

(c)  Separating  approximately  equal  numbers  of  points* 

The  slopes  and  heights  to  -which  these  lines  may  be  considered  representa- 
tive are  shovm  in  Table  11  for  each  sounding* 

laximum  height     Values  of  K]_ 
Sounding                       of  Ij-near  distri-   in  equation  (5  •5) 
(Time  in  LCT)   Slope  of  line,  m   bution  (feet) (cm^"^  seo"*^) 

0200  1.14  23.5  565 

0500  1.17  23.5  506 

1100  1.79  45.0  556 

1300  1.69  45.0  2660 

1600  1,70  23.5  410 

1900  1.05  45.0  155 

2200  1,23  16,0  M4 

Slopes  of  log  K  Against  Log  Z  for  the  Lines  of  Plate  I 


TABLE  11. 
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There  is  a  tendency  for  the  lines  for  the  1^00  and  2200  LCT  soiindings 
to  continue  to  be  representatire  up  to  2^5  feet,  Ho'vTever,  for  the  0200  and 
0500  LCT  soundings  beyond  25«5  feet,  the  points  fall  well  belo--;/  the  line. 
This  suggests  the  possibility  of  fitting  a  curve  other  than  a  straight  line 
(perhaps  a  parabola)  in  the  early  morning  hours;  hov/ever,  this  idoa  v/as  not 
pursued  further. 

If  vre  vnrite  the  equation  of  any  straight  line  of  Plate  I  in  the  form 

log  K  -  m(log  Z)  +b,  (5.2) 

■where  m  is  the  measured  slope    (see  Table  11)  of  the  line,  vre  obtain  upon 
raising  to  povrers  of  10,  the  equation 

K  5  10^  Z^  s  Kj  Z^,  (S3) 

Since  each  line  passes  through  the  point  corresponding  to  Z  s  I/4  foot,  the 
value  of  K^  for  each  sounding  may  be  assigned  by  setting  Z  -  I/4  and  K  to 
its  appropriate  value.  The  values  of  K-j_  for  each  sounding  are  shoTm  in  the 
last  column  of  Table  11.  Equation  (5 •J)*  with  the  appropriate  values  of  Kj_ 
suid  m,  represent  the  distribution  of  K  vfith  respect  to  height  for  the  various 
times  of  day  and  up  to  the  heights  indicated  in  Table  11.  This  is  the  chief 
empirical  result  of  this  study. 

Finally,  -we  shall  state  the  conclusions  resulting  from  the  foregoing 
analysis  and  from  Figures  4  through  8.  These  conclusions  are: 

(1)  The  average  mid-day  value  of  m  is  1.75* 

(2)  The  average  value  of  m  near  midnight  is  l.l8. 

(5)  The  value  of  m  seems  to  approach  unity  at  sunset. 
(4)  The  exponent  m  appears  to  be  a  function  of  stability. 
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(5)  The  inaxiinum  value  of  K,  level  for  level  up  to  60  feet,  appears 
to  occur  at  about  1^00  LCT,  at  about  the  time  of  the  maximum  surface  tem- 
perature. The  povrer  law  satisfied  by  K  at  this  time  of  day  is 

K  z   2660  Z^*^^, 
where  Z  is  the  height  in  feet» 

(6)  The  minimum  value  for  K,  level  for  level  up  to  60  feet,  appears 

to  occur  about  I9OO  LCT,  with  a  very  slcfw  rate  of  increase  during  the  night* 

(7)  From  Figure  8  at  125  feet  there  appears  to  be  a  second  pronounced 
diurnal  maximj.xm  of  K  at  2200  LCx#  This  appears  to  be  due  to  the  diurnal 
variation  in  the  height  of  the  surface  layer*  According  to  Lettau  [Sj  « 
the  surface  la^,'er  is  characterized  by  increasing  values  of  K  v/ith  height, 
with  decreasing  values  of  K  occurring  above  this  layer*  If,  during  the  day- 
li£;ht  hours  the  125  foot  level  is  above  the  surface  layer,  but  passes  into 
it 'at  about  2200  LCT,  the  large  increase  in  the  value  of  K  at  125  feet  ap- 
pears to  be  reasonable* 
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